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1. 

FNWC 

Fleet  Numerical  Weather  Central,  Monterey,  California  (now  Fleet 
Numerical  Oceanography  Center) 

2. 

IGDR 

- 

Interim  Geophysical  Data  Records 

3. 

SDR 

- 

Sensor  Data  Record 

4. 

PGD 

- 

Filtered  Geophysical  Data 

5. 

UGD 

- 

Unfiltered  Geophysical  Data 

6. 

SGD 

- 

Segmented  Geophysical  Data 

7. 

NASA 

- 

National  Aeronautics  and  Space  Administration 

8. 

PMDF 

- 

Project  Master  Data  File 

9. 

rev 

- 

One  revolution  of  the  satellite  around  the  earth  commencing  at  the 
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INTRODUCTION 


SEASAT,  a  follow-on  of  GEOS-3,  was  a  National  Aeronautics  and  Space 
Administration  (NASA)  satellite  for  measuring  global  ocean  dynamics  from  space. 
The  spacecraft  was  launched  2t>  June  1978  with  instruments  aboard  that  provided 
data  on  wave  height,  wind  speed  and  direction,  ice  fields,  ocean  surface  topog¬ 
raphy  and  weather.  Data  over  the  oceans  was  stored  aboard  and  transmitted  when 
in  the  range  of  the  observation  stations  located  at  Merritt  Island,  Florida, 
Madrid,  Spain,  and  Fairbanks,  Alaska.  Goddard  Space  Flight  Center  in  Greenbelt, 
Maryland  obtained  the  data  collected  by  the  3  stations  and  distributed  it  to 
the  Jet  Propulsion  Laboratory  (JPL)  in  Pasadena,  California,  and  Fleet 
Numerical  Weather  Center  (FNWC)  in  Monterey,  California,  where  data  was  con¬ 
verted  to  engineering  units  and  some  preprocessing  done.  Early  in  the  mission, 
data  from  Fairbanks  was  transmitted  by  telecommunications  satellite  to  FNWC, 
but  this  process  was  discontinued  because  of  an  operational  problem.  The  mis¬ 
sion  ended  9  October  1978  with  approximately  1000  revolutions  (revs)  of  data 
collected.  As  mentioned  earlier,  the  spacecraft  carried  a  number  of  different 
instruments;  however,  this  report  is  only  concerned  with  the  radar  altimeter. 

The  purpose  of  the  radar  altimeter  was  to  measure  significant  wave  height 
at  the  subsatellite  point  and  the  precise  altitude  between  the  satellite  and 
the  ocean  surface.  The  altitude  when  combined  with  the  Doppler  precise  orbit 
(Colquitt  et  al.,  1980],  atmospheric  and  environmental  effects,  gives  the  sea- 
surface  topography,  which  is  an  approximation  of  the  geoid  height.  This  report 
describes  the  SEASAT  altimeter  data-processing  system  which  employs  the  alti¬ 
tude  measurements  in  the  determination  of  best-estimated  along-track  geoid 
heights  and  vertical  deflections. 


SEASAT  ALTIMETER  DATA-PROCESSING  SYSTEM 


GENERAL  DESCRIPTION 

The  Naval  Surface  Weapons  Center  (NSWC)  altimeter  data-processing  system 
entails  the  following:  the  generation  of  the  Interim  Geophysical  Data  Record 
( IGDR)  tapes;  separation  of  revs  and  time-continuous  segments  of  revs;  computa¬ 
tion  of  sea-surface  heights;  and  application  of  an  adaptive  Kalman  smoother  to 
obtain  best-estimated  along-track  geoid  heights  and  vertical  deflections.  An 
overview  of  the  processing  system  is  shown  in  Figure  1.  At  FNWC,  Project 
Master  Data  File  (PMDF)  tapes  are  used  in  the  generation  of  Sensor  Data  Record 
(SDR)  tapes,  which  are  then  used  as  input  to  NSWC  software  (REDUCE).  REDUCE 
outputs  IGDR  tapes  which  were  shipped  from  FNWC  to  Dahlgren.  Formats  for  the 
IGDR  and  SDR  tapes  and  the  REDUCE  program  are  described  in  an  intralaboratory 
Technical  Note  (80-96)  by  J.  M.  Futcher,  Jr.  In  brief,  the  REDUCE  program 
separates  land  and  ocean  data,  computes  atmospheric  and  environmental  correc¬ 
tions,  and  prefilters  and  aggregates  the  data.  The  altimeter  data  given  on  the 
SDR  at  10  points  per  second  (pps)  can  be  aggregated  by  REDUCE  to  a  rate  speci¬ 
fied  by  an  input  parameter.  Data  was  aggregated  to  2  pps  for  the  overall 
SEASAT  altimeter  data  reduction. 
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Figure  1.  SEASAT  Radar  Altimeter  Data-Processing  System 


The  IGDR  tapes  are  the  major  input  to  the  altimeter  data-processing  at 
NSWC.  The  tapes  contain  such  information  as  rev  number,  time,  latitude,  longi¬ 
tude,  altimeter  range,  and  atmospheric  and  environmental  corrections.  Each 
tape  contains  1  to  8  data  files  where  each  file  corresponds  to  a  dump  of  the 
recorder  on  board  the  satellite.  A  file  is  made  up  of  records,  one  for  each 
data  point.  Data  within  a  satellite  rev,  revs  within  a  file,  and  most  files 
are  time  ordered. 

Data  files  created  during  the  processing  are  the  Log  Index,  Segmented 
Geophysical  Data  (SGD),  Unfiltered  Geophysical  Data  (UGD) ,  and  Filtered 
Geophysical  Data  ( FGD) .  The  SGD  and  UGD  files  are  temporary  but  can  be  saved 
for  special  purposes  through  program  control. 

The  FGD  file  is  the  major  output  from  the  altimeter  data-processing  system 
and  contains  basically  sea-surface  heights,  corrections,  best-estimated  along- 
track  geoid  heights  and  vertical  deflections.  Both  the  IGDR  and  FGD  are  main¬ 
tained  at  NSWC  and  have  been  distributed  to  other  specified  agencies  for  analy¬ 
ses  and  for  determining  other  geodetic  parameters.  The  next  4  sections  will 
contain  brief  descriptions  of  the  major  programs  used  in  the  processing. 


SEGEMENTATION  OF  DATA  INTO  REVS 


The  function  of  the  Segmenter  program  entails  separation  of  data  from  the 
IGDR  tapes  into  revs,  division  of  revs  into  time-continuous  segments,  calcula¬ 
tion  of  tide  corrections,  determination  of  some  environmental  corrections,  and 
creation  of  the  Log  Index  and  the  Segmented  Geophysical  Data  files.  Data  for 
each  rev  is  divided  into  time-continuous  segments  and  processed  as  a  separate 
entity.  Time-continuous  segments  are  determined  by  examining  the  rev  data  for 
the  existence  of  a  constant  time  interval  (At  -  .5  sec) .  If  the  time  interval 
is  not  approximately  .5  sec,  but  <5  sec,  a  linear  interpolation  procedure  is 
used  to  dub  in  the  number  of  points  required  to  bridge  the  interval.  If  the 
interval  is  >5  sec,  the  previous  time-continuous  segment  ends  and  a  new  segment 
begins. 
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For  each  altimeter  measurement  the  ocean  tide  correction  is  determined  by 
the  Schwiderski  Ocean  Tide  model  (1979).  This  model  consists  of  the  7  leading 
tide  components:  M2 ,  S2 ,  N2 ,  K2 ,  K^ ,  Oj ,  and  Pj  of  the  solid  earth,  oceans,  and 
the  tidal  yielding  of  the  solid  earth.  Such  information  as  the  time  of  initial 
point  in  segment,  time  interval  between  points,  longitude  of  the  ascending 
node,  equator  crossing  time,  and  number  of  points  in  a  segment  is  required  by 
the  model. 

The  timing  correction  (Appendix  B)  is  applied  to  the  observation  time  in 
the  Segmenter  and  corrections  for  the  tilt/sea  state,  geometric  factor,  elec¬ 
trical  delay  distance  are  computed.  Algorithms  for  these  corrections  are  given 
in  Appendix  B. 

In  the  altimeter  data,  erroneous  values  of  latitude  and  longitude  were 
found  periodically.  When  these  values  were  encountered  corresponding  correc¬ 
tions  dependent  upon  latitude  and  longitude  were  set  to  zero  or  a  standard 
value. 


The  Log  Index  file  created  by  the  segmenter  is  used  as  input  to  an  index 
program,  which  provides  descriptive  information  for  each  rev  and  segment  of  rev 
in  the  complete  SEASAT  data  base.  The  SGD  file  is  used  as  input  to  the  second 
major  program,  GENERATE,  which  determines  the  sea-surface  heights. 


COMPUTATION  OF  SEA-SURFACE  HEIGHTS 


The  height  of  the  a ltimetric  geoid  relative  to  the  reference  ellipsoid 
(sea-surface  height)  is  obtained  by  combining  the  altimeter  range  with  the 
Doppler  precise  orbit  and  applying  environmental  and  atmospheric  effects.  The 
time  of  the  altimeter  range  is  used  in  an  8-point  Lagrangian  interpolation  to 
compute  subsatellite  position,  latitude  and  longitude.  Algorithms  for  deter¬ 
mining  the  subsatellite  position  and  sea-surface  height  were  provided  by  Ray 
Manrique  and  are  given  in  Appendix  A.  Corrections  applied  to  the  altimeter 
range  are  shown  in  Figure  2. 

The  computer  program,  GENERATE,  performs  the  computation  of  the  sea-sur¬ 
face  height  and  outputs  the  temporary  UGD  file  for  use  in  the  filtering  section 
of  the  data-processing  system. 


GENERATION  OF  BEST-ESTIMATED  ALONG-TRACK  GEOID  HEIGHTS  AND 
VERTICAL  DEFLECTIONS 


The  sea-surface  heights  are  reduced  to  best-estimated  along-track  geoid 
heights  and  vertical  deflections  by  a  Kalman  smoother  based  upon  the  third- 
order  Markov  process  (Jordan,  1972).  The  Kalman  smoother  for  SEASAT  was  pro¬ 
vided  by  Dr.  C.  J.  Cohen  and  is  described  in  West  et  al.  (1977).  An  adaptive 


3 


CENTER  GRAVITY  SEASAT 
Oq.  GEOMETRIC  DISTANCE 
ANTENNA  ELECTRICAL  CENTER 


HEIGHT 
ABOVE 
REFERENCE  « 
ELLIPSOID, 
h 


> 


TRUE  RANGE  RADAR  MEASUREMENT  |H| 

DRY  TROPOSHERE  DISTANCE  EQUIVALENT.  Dp 
WET  TROPOSHERE  DISTANCE  EQUIVALENT.  Dw 
IONOSPHERE  DISTANCE  EQUIVALENT.  D, 

RADAR  INSTRUMENT  DELAY  DISTANCE 
EQUIVALENT.  DE 

TILT/SEA  STATE  DISTANCE  EQUIVALENT.  DTSWH 


ELLIPSOID 


GH  h  H  D6  ♦  De  ♦  Dd  *  Ow  ♦  0, 
‘  °T  '  °B  *  °TSWH  '  DS 


1.  HIGH  TIDE  (Dt)  ARE  ASSUMED  POSITIVE  (*) 

2  BAROTROPIC  (DB)  ARE  POSITIVE  (*). 

IF  P  ■  1013  mb 

3  CURRENT  CORRECTION  OMITTED 

4  SEA  STATE  CORRECTION.  DswH  ,s  COMBINED  WITH 


Figure  2.  Atmospheric  and  Environmental  Corrections 


TSWH 


4 


algorithm  written  by  Ugincius  (1977)  is  useii  +01  improvement  of  a  priori  para¬ 
meters;  auto-correlation  distance,  guord-height  variance,  and  the  noise  stan¬ 
dard  deviation  required  by  the  smoother.  The  autocorrelation  distance  is  de¬ 
termined  for  ’Dull  km  spans  of  data  and  spans  with  similar  autocorrelation 
values  are  grouped  as  a  segment. 

Magnitudes  of  the  geoid  heights  and  vertical  deflections  sometimes  fall 
outside  the  expected  range.  These  values  may  be  caused  by  such  factors  as  loss 
of  bits,  and  the  presence  of  land  data  near  the  coastlines  anti  over  small  land 
masses.  Geoid-height  values  less  than  -125  m  and  greater  than  +125  m  are  set 
to  -125  m  or  +125  m,  respectively.  Vertical  deflections  are  handled  similarly, 
those  values  less  than  -100  arc  sec  and  greater  than  +100  arc  set  are  set  to 
-100  arc  sec  or  +100  arc  sec,  respectively.  If  either  of  the  two  conditions 
exist  the  longitude  value  is  made  negative  to  indicate  that  data  has  been 
dubbed  in.  A  flag  word  in  the  FGD  file  is  also  set  to  reflect  this  data.  In 
the  smoothing  procedure,  dubbed-in  data  points  have  zero  weights. 

Output  from  the  Kalman  smoother  program  is  the  FGD  file,  which  consists  of 
filtered  data  from  one  satellite  rev.  Such  information  as  sea-surface  heights, 
corrections  applied,  statistical  information,  and  best-estimated  along-track 
geoid  heights  and  vertical  deflections  make  up  the  FGD  file.  The  rev  data  base 
consists  of  over  900  FGD  files  on  20  device  sets.  Each  FGD  file  is  a  permanent 
file  with  file  names  of  the  form  FGDRRRRR  where  RRRRR  is  the  rev  number.  The 
rev  data  base  was  used  in  the  generation  of  an  area  data  base  that  consists  of 
48  permanent  files  on  9  device  sets.  The  area  base  differs  from  the  rev  base 
in  that  erroneous  data  points  at  both  ends  of  a  segment  are  eliminated  and  the 
data  has  been  aggregated  to  1  pps  and  data  in  an  area  only  contains  the  portion 
of  the  rev  that  passes  through  the  area.  The  surface  of  the  earth  is  divided 
into  areas  as  shown  in  Figure  3.  Area  file  names  are  of  the  form  FGDSSAREAXX 
where  SS  means  SEASAT  and  XX  is  the  area  number. 


PLOT  OF  FILTERED  GEOPHYSICAL  DATA  FILE 


Parameters  from  the  FGD  file  may  be  plotted  on  the  SD4060  or  Tektronix 
(4051)  plotting  equipment.  Parameters  which  can  be  plotted  versus  time  are  the 
following:  sea-surface  heights,  geoid  heights  (Figure  4),  vertical  deflections 

(Figure  5),  significant  wave  heights  (Figure  6),  and  the  automatic  gain  control 
(Figure  7).  A  complete  satellite  rev  may  be  shown  on  one  plot  frame  or  divided 
into  time  spans  and  each  time  span  shown  on  one  frame. 


CONTROL  PROGRAM 

Reduction  of  the  altimeter  data  is  controlled  by  the  SEASAT  Radar 
Altimeter  Data-Processing  System  (SRAPS)  computer  program  (J.  D.  Clark,  unpub¬ 
lished  user's  guide).  SRAPS  gives  the  user  the  flexibility  in  a  single  comput¬ 
er  job  to  control  the  execution  of  the  set  of  computer  programs  required  to 
reduce  the  SEASAT  altimeter  measurements,  taken  over  t’.ie  oceans,  to  filtered 
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along-track  geoid  heights  and  vertical  deflections.  The  computer  programs  are 
executed  in  the  order  described  in  Figure  1.  A  file  containing  the  information 
for  reducing  the  altimeter  data  is  given  to  the  SRAPS  program,  which  generates 
job  control  necessary  for  processing  a  set  of  data.  The  job  control  is  then 


Figure  3.  Filtered-Area  Data  Base 


Figure  7.  Automatic  Gain  Control  for  Rev  1348 
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Figure  A-1.  Height  of  Geoid  Above  Reference  Ellipsoid 


h  -  altimetric  height  above  the  geoid  (km) 
rs  -  distance  of  satellite  from  center  of  reference  ellipsoid  (km) 
N  -  height  of  geoid  above  reference  ellipsoid  (m) 

<J>  -  geodetic  latitude  (deg) 

H  -  geometric  height  above  the  ellipsoid  (km) 
a  -  semi-major  axis  of  reference  ellipsoid  (km) 
e  -  eccentricity  of  the  reference  ellipsoid 
N  =  H  -  h 


where 


a( 1-e  ) 


sin4> 


e2  sin2$ 


and 


h  is  corrected  for  atmospheric  and  environmental  effects  as  shown  earlier  in 
Figure  2. 


4> i  =  <J>k+1  and  X,  are  computed  for  each  observation  time,  t^. 
a  =  6378.1^*5  km,  e2  =  2f  -  f2,  and  f  =  1/298.25 
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APPENDIX  B 


ATMOSPHERIC  AND  ENVIRONMENTAL  CORRECTIONS 


1.  Tilt/ Sea-state  Correction  to  Height 

The  Wallops  Flight  Center  Tilt/ Sea  state  Correction  Tables  (B.  F. 
Townsend,  3/26/79)  were  fitted  with  the  following  bi-linear  (Anderle 
function) : 

AH(cm)  =  -5.  +  9.33{T)  +  1.0625(SWH)  +  1 .75(T) ( SWH) , 
where:  T  =  tilt  in  degrees, 

SWH  =  significant  wave  height  in  meters. 

The  maximum  deviation  from  the  tabular  values  is  ~  10  cm. 

The  correction  is  labeled  in  the  printed  output  as  Tilt/ SWH  Correction  to  H. 

2.  The  SWH  correction  obtained  from  a  similar  fit  is  given  below. 

1/3  2 

AH  (m)  =  A SWH  =  -.32  x  SWH  x  T2 

where,  T  =  tilt  in  degrees  and  AH1//3  is  given  in  meters. 

3.  Wet  Tropospheric  Correction 

The  wet  trospheric  correction,  Djj,  is  computed  from  input  surface 
temperature  and  partial  pressure  of  water  vapor  given  below. 

Dw(cm)  =  86400.  Pw/(273.  +  Tg)2 , 

where 

Pw  =  partial  pressure  of  water  vapor  (mB) , 

Tg  =  surface  temperature,  deg.  C. 

If  Pjj  is  unavailable,  use  12.272. 

If  Tg  is  unavailable,  use  0. 

4.  Dry  Tropospheric  Correction 

The  dry  tropospheric  correction  Dp,  is  computed  from  input  pressure  and 
latitude  as  given  below. 

Dn(cm)  =  <2-277  -  P, 

D  10 
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where 


$  =  latitude  from  location  data, 

P  =  pressure  (mB). 

If  P  is  unavailable,  use  1013.3. 

5.  Barotropic  Correction 

The  barotropic  correction,  Dg,  computed  from  input  atmospheric  pressure 
as  given  below. 

DB(m)  =  -. 009948 (Pa  -  1013.0), 

where 

pa  *  atmosphereic  pressure  (mB)  supplied  by  FNWC. 

6.  Ionospheric  Refraction  Correction 

The  ionospheric  refraction  model  shown  in  Figure  B-l  (Bent,  et  al.,  1975) 
was  developed  by  James  Clynch  and  Arnold  Tucker  of  Applied  Research  Laboratory, 
University  of  Texas,  Austin,  Texas  and  modified  by  Dr.  Ralph  Gibson  of  NSWC. 

The  model  uses  telecommunications  predictions  obtained  from  the  Environmental 
Data  Service  (now  known  as  National  Oceanic  and  Atmospheric  Administration). 

The  ionospheric  correction,  AR,  is  given  by  the  equation  (1)  below. 


Figure  B-l.  Ionospheric  Refraction  Model 

For  SEASAT  altimeter  k=0,  since  the  range  is  straight  up  and  down;  therefore 
the  ionospheric  refraction  correction,  Dx ,  is  given  by  the  equation  (2)  and 
ranges  between  3  and  15  cm  with  a  ±  15%  accuracy. 


20 


r 


D  (cm)  =»  —  /  N(x,y,z)  dr, 

tL  r 

e 

where 

r  =  distance  from  center  of  earth  to  the  satellite  (km) 
s 

r0  =  the  radius  of  the  earth  (km) 
f  =  transmitted  frequency  (MHz) 

C  *  a  known  constant  (=  8.061389  x  10“5) 

N(x,y,z)  =  electron  density  (el/cm3) 

The  vertical  electron-density  profile  is  shown  in  Figure  B-3.  The  E  and  Fl 
layers  are  modeled  by  a  parabola;  F2  layer  up  to  1^2  is  modeled  by  a  fourth 
degree  curve;  H„,F2  to  Hcrit  is  modeled  by  a  parabola;  and  the  remaining 
3  portions  of  F2  are  modeled  by  a  decaying  exponential  with  decay  con¬ 
stants  Kj ,  K2 ,  and  Kj . 

PARAMETERS  FOR  IONOSPHERIC  ELECTRON -DENSITY  PROFILE 

-  Critical  frequency  of  E  region 

-  Height  of  maximum  ionization  of  E  region 

-  Semithickness  of  E  region 

-  Critical  frequency  of  Fl  region 

-  Height  of  maximum  ionization  of  Fl  region 

-  Semithickness  of  Fl  region 

-  Critical  frequency  of  F2  region 

-  Height  of  maximum  ionization  of  F2  region 

-  Semithickness  of  lower  and  middle  F2  region 
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o 


H  F2 
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K.j  ,K^  »Kj  -  Decay  constants  for  upper  F2  region  (topside) 


Hcr^t  -  Height  at  which  middle  and  upper  F2  regions  match  gradients 

H1m  -  Center  of  lower  topside 

H2m  -  Center  of  Middle  Topside 

H3 

m  -  Center  of  Upper  Topside 

7.  Geometric  and  Radar  Instrument  Delay  Distance 


a)  The  sum  (d  +  dE)  given  for  the  geometric  factor  has  the  electrical 
delay  built  in.  For  SEASAT  altimeter  processing,  geometric  distance,  D^,  is 
set  to  d  +  dg  and  radar  instrument  delay  distance,  D£,  is  set  to  zero.  The 
algorithm  for  d  as  shown  below  was  obtained  from  H.  Hagar  (unpublished  data, 
1978)  at  the  Jet  Propulsion  Laboratory,  California  Institute  of  Technology, 
Pasadena,  California. 


d(m)  =  .0254 


'1086187 


434.5  W. 


4852.6 


■33.6 


CENTER  OF  GRAVITY 


d 


ALTIMETER 

BASEPLATE 


ALTIMETER 

ELECTRICAL 

CENTER 


P’igure  B-4.  Geometric  and  Radar 
Instrument  Delay  Distance 
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where 


Wh  =  104.5  -  WhR  (lbs) 
w 

hR  =  pounds  of  hydrazine  remaining 
WhR  =  89.4  lbs  on  13  July  1978 

Wj,R  is  updated  after  each  major  orbit  change/momentum  wheel  dump.  These 
values  are  given  in  Table  B-1. 

Table  B-1.  Pounds  of  Hydrazine  Remaining,  W^R,  Year  1978 
Maneuver 


nuuc u» 

Date 

Day 

WhR 

6/27 

178 

89.4 

8/15 

227 

87.97 

8/18 

230 

86.05 

8/23 

235 

84.69 

8/26 

238 

75.17 

9/10 

252 

74.62 

b)  dE(m)  =  0.47752  +  4.57x10-9  x  c, 

where 

c  =  speed  of  light  in  meters/second, 
c  =  299,792,458  m/s  (1974  value) 

dg  was  obtained  from  J.  Lorell  (unpublished  data,  1978)  at  the  Jet  Propulsion 
Laboratory. 

8.  Ocean-Tide  Correction 

The  ocean-tide  correction  was  computed  by  the  Schwiderski  ocean-tide  model 
(1979),  which  comprises  the  seven  (7)  leading  tide  constituents;  M2 ,  S2 ,  N2 , 

K2 ,  ,  Oj ,  and  Pj^  of  the  solid  earth,  oceans,  and  the  tidal  yielding  of  the 

solid  earth.  The  total  instantaneous  ocean  tide  is  computed  with  10  cm  accura¬ 
cy.  The  tide  reduction  program  is  based  on  a  one-by-one-degree  amplitude  and 
phase  tide  table,  which  was  hydrodynamics lly  computed  by  an  improved  version  of 
the  numerical  model  by  W.  Zahel  (1977). 


9.  Timing  Correction 


In  an  orbit  accuracy  assessment  for  SEASAT  by  Schutz  and  Tapley  (1980)  an 
altimeter  time-tag  correction  of  -.079  sec  was  recommended.  This  value  differs 
by  -.028  sec  from  the  value  -.051  sec  derived  earlier  from  intersection  analy¬ 
sis  of  SEASAT  data  at  NSWC/DL.  The  time-tag  correction  -.051  sec  was  applied 
to  the  altimeter  data  during  reduction  to  along-track  geoid  heights  and  verti¬ 
cal  deflections.  A  time-tag  correction  was  not  applied  to  the  altimeter  data 
at  FNWC. 
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